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ABSTRACT 

We have studied the X-ray luminosity function (XLF) of low-mass X-ray binaries (LMXBs) in 
the nearby lenticular galaxy NGC 3115, using the Megasecond Chandra X-Ray Visionary Project 
Observation. With a total exposure time of ~1.1 Ms, we constructed the XLF down to a limiting 
luminosity of ~10 36 erg s~ 4 , much deeper than typically reached for other early-type galaxies. We 
found significant flattening of the overall LMXB XLF from dN/dL oc L _2 ’ 2± °' 4 above 5.5 x 10 37 
erg s _1 to dN/dL oc £- 10 ± 01 below it, though we could not rule out a fit with a higher break at 
~1.6 x 10 38 erg s -1 . We also found evidence that the XLF of LMXBs in globular clusters (GCs) is 
overall flatter than that of field LMXBs. Thus our results for this galaxy do not support the idea 
that all LMXBs are formed in GCs. The XLF of field LMXBs seems to show spatial variation, with 
the XLF in the inner region of the galaxy being flatter than that in the outer region, probably due 
to contamination of LMXBs from undetected and/or disrupted GCs in the inner region. The XLF 
in the outer region is probably the XLF of primordial field LMXBs, exhibiting dN/dL oc L _1 ‘ 2±0-1 
up to a break close to the Eddington limit of neutron star LMXBs (~1.7 x 10 38 erg s _1 ). The break 
of the GC LMXB XLF is lower, at ~1.1 x 10 3 ‘ erg s _1 . We also confirm previous findings that the 
metal-rich/red GCs are more likely to host LMXBs than the metal-poor/blue GCs, which is more 
significant for more luminous LMXBs, and that more massive GCs are more likely to host LMXBs. 
Subject headinqs: X-rays: binaries — globular clusters: general - Galaxy:stellar content — X-rays: 
individual (NGC 3115) 


1. INTRODUCTION 

Population studies of X-ray binaries in nearby galax¬ 
ies have been made possible thanks to the superb spa¬ 
tial resolution and excellent sensitivity of the Chandra 
X-ray Observatory (|Weisskopf et all 12002 ). The X-ray 
luminosity functions (XLFs) of point sources have been 
obtained for many galaxies, and they are found to be 
environment dependent (see iFabbianol 120061 for a re¬ 
view). In young normal galaxies, high-mass X-ray bi¬ 
naries (HMXBs) dominate, and the XLFs follow a sim¬ 
ple power law (PL) dN/dLx oc L/j/* with a Ri 1.6 over 
a large r ange of luminosity: 10 35 erg s~ 4 < Lx < 10 40 
erg s _1 dGrimm et all 12003 : iMineo et al.l I2012T) . In old 
normal galaxies or in the bulge of young normal galax¬ 
ies, low-mass X-ray binaries (LMXBs) dominate, and the 
XLFs seem relatively complicated, showing bot h a high- 
luminosity break at Lx ~ 5 x 10 38 erg s~ 1 2 3 4 (ISarazin et ah 
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l2012f) and a low-luminosity break at Lx ~ 5 x 10 37 erg 
s~ 4 (iGilfanovI 12004 IVoss et all 120091 IKim et all 120091 : 
IZhang et alJl2012t ILuo et al.l 120121) "Die slope between 
these two breaks is a rU.8-2.2. Above the high- 
luminosity break, the XLFs decrease sharply. Below the 
low-luminosity break, the XLFs might flatten to a « 1.0. 

The high-luminosity break may be due to the Ed¬ 
dington limit of neutron star (NS) LMXBs. The 
low-luminosity break has been attributed to either 
the transition from mass transfer driven by magne¬ 
tized stellar wind at high luminosities to mass trans¬ 
fer driven by gravitat ional wave e mission at low lumi¬ 
nosities dPostnov fe Kuranovll2005ll . or different types of 
donor stars with the high-luminosity ones being giants 
and the low-luminosity ones being main-sequence stars 
(jRevnivtsev et all 1201 lfl . The former predicts the slope 
below the break to be about 1.0 and the slope above the 
break to be about 2.0. For the latter explanation, the 
steepening of the XLF at high luminosity is due to the 
short life time of binary systems with giants. 

The XLFs of LMXBs have often been obtained by com¬ 
bining multiple galaxies in order to improve the statis¬ 
tics. However, this method is subject to the limitation 
that the normalizations of the XLFs in diffe rent g alaxies 
show a scatter of more than a factor of two (IZhang et al.l 
(2012 ). The XLFs of LMXBs well below 10 37 erg s _1 are 
still only obtained for very fe w old galaxies, most notably 
Centaurus A and NGC 3379 (IVoss et al.ll2009t [Kim et al.l 
I2009H . and the bulge of M31 ( Voss fe Gilfanovll2007allbH . 
A low detection limit is critical for constraining the low- 
luminosity break, which has been shown m ostl y for the 
old populations in young galaxies (lFabbianoll2006H . Deep 











































2 


Lin et al. 



Fig. 1.— Chandra X-ray image of NGC 3115 in 0.5-7 keV. The 
image is adaptively smoothed with the CIAO task csmooth and 
exposure corrected. The D 25 ellipse of the galaxy and the approx¬ 
imate FOV of the HST /ACS mosaic observation (dashed box) are 
also shown. 

observations of old galaxies are also needed for the inves¬ 
tigation of the differences between the XLFs of LMXBs 
in globular clusters (GCs) and in the stellar field, which 
can be used to check whether they have the same ori¬ 
gin. For instance, there is a relative underabundance of 
faint LMXBs in GCs when compared with field LMX Bs 
(jVoss et al.ll2009t iKim et all 1 2 0091 iZhang et all 1 2 01 il l. 

NGC 3115 was selected as the target of a 1 Megasec- 
ond Chandra X-ray Visionary Project (XVP) in Cycle 13. 
One main goal was to study the gas flow inside the Bondi 
radius of the central supermassive black hole (BH), which 
has been reported in lWong et all ( 2014 ). The other goal 
was to have a deep look at the X-ray binary population 
of a normal early-type galaxy. The detailed analysis of 
the data, including the source list and detail ed properties 
of special sources, will be presented in iLin et all (12015L 
Paper I hereafter). In the present paper, we concentrate 
on the XLF of LMXBs, especially on its faint end below 
10 37 erg s _1 . NGC 3115 is a lenticul ar ( SO) galaxy with 
an age of 8.4 ±1.1 Gyr (iSancliez-Blazciuez et al.l 120068) 
and at a distance of 9.7 Mpc (|Tonrv et al.1 120011) . In¬ 
cluding previous observations, the total exposure time of 
Chandra on this galaxy is ~1.1 Ms, reaching a limiting 
luminosity of ~ 10 36 erg s _1 . Thus it is one of the best 
observed normal early-type galaxies by Chandra. Ac¬ 
companying the Chandra XVP observation, there is a 
six pointing Hubble Space Telescope (HST) mosaic ob¬ 
servation in the F475W and F850LP filters (hereafter g 
and z filters, respectively) using the Advanced Camera 
for Surveys (ACS). It provides th e information of GCs in 
the galaxy (IJennings et alJl20l4 h. which will be used by 
us to investigate the dependence of the XLF of LMXBs 
on the stellar environment. 

The paper is organized as follows. In Section [2] we 
describe the X-ray data reduction, construction of XLFs, 
incompleteness correction, and cross-correlation of the 


TABLE 1 
Observation Log 


Notation 

Obs. ID 

Date 

Exposure 

(ks) 

Offset a 

(arcmin) 

1 

2040 

2001-06-14 

35.8 

1.5 

2 

11268 

2010-01-27 

40.6 

0.1 

3 

12095 

2010-01-29 

75.6 

0.1 

4 

13817 

2012-01-18 

171.9 

0.0 

5 

13822 

2012-01-21 

156.6 

0.0 

6 

13819 

2012-01-26 

72.9 

0.0 

7 

13820 

2012-01-31 

184.1 

0.0 

8 

13821 

2012-02-03 

157.9 

0.0 

9 

14383 

2012-04-04 

119.4 

0.3 

10 

14419 

2012-04-05 

46.3 

0.3 

11 

14384 

2012-04-06 

69.7 

0.3 


a Aim point offset from observation 13820. 


X-ray and optical sources. In Section [3l we show the 
spatial distributions of different populations, present the 
total LMXB XLF, compare the XLFs of GC and field 
LMXBs, and investigate the GC LMXB properties. In 
Section^ we discuss various caveats on the XLFs that we 
obtain and the implication of our results for the nature 
and formation of LMXBs. Our conclusions are given in 
Section [5] 

2. DATA ANALYSIS 

2.1. Observations and Source Detection 

The Chandra observations of NGC 3115 are listed in 
Table ED There are 11 observations in total from essen¬ 
tially three epochs: one in 2001, two in 2010 and nine 
in 2012. All observations used the imaging array o f the 
AXA F CCD Imaging Spectrometer (ACIS; iBautz et a. 1.1 
119981) . The reduction of the data and the creation of the 
source list were presented in detail in Paper I, and here 
we briefly summarize the procedure adopted. The data 
were analyzed with the Chandra Interactive Analysis of 
Observations (CIAO, version 4.6) package. The data 
were reprocessed to apply the latest cali brat ion jCALDB 
4.5.9) and the subpixel algorithm dLi et al.l 120041) using 
the CIAO script chandra-repro. Some short background 
flares seen in observations 2040, 13819 and 13822 were 
excluded. The final exposure used for each observa¬ 
tion is given in Table [T] The relative astrometry be¬ 
tween observations was corrected, and the source detec¬ 
tion was performed on individual observations as well 
as on the merged one using the 0.5-7 keV energy band 
with the C IAO wavdetect wavele t-based source detection 
algorithm dFreeman et al.l [20 02b We used two different 
image binning resolutions: one at single sky pixel reso¬ 
lution (O'/492) over the full field of view (FOV) and the 
other at 1/8 sky pixel resolution covering an area of 3'x3' 
centered at the center of NGC 3115. The subpixel bin¬ 
ning images were used to improve the spatial resolution 
of the crowded field near the center of the galaxy. Sources 
detected from the merged observation and individual ob¬ 
servations were cross-correlated to create the final unique 
source list. 

For each unique source, we extracted the source and 
background spectra and created the response file for each 
individual observation. They were then merged to cre¬ 
ate the spectra and response files for the merged obser¬ 
vation. The source region was set to be a circle enclos¬ 
ing 90% of the point spread function (PSF) at 2.3 keV. 
The background region was set to be a concentric annu- 
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lus, with inner and outer radii of two and five times the 
source radius, respectively. The background-subtracted 
0.5-7 keV count rates were then converted to unabsorbed 
fluxes and luminosities, with the conversion factors based 
on the corresponding response files and assuming an ab¬ 
sorbed PL spectral shape with a photon index of 1.7 
and the Galactic ab sorption IVh = 4.32 x 10 20 cm ~ 2 
dKalberla et al. I [20051) . 

2 . 2 . Incompleteness Calculation and XLF Construction 

The point-source detection sensitivity varies across the 
Chandra image, owing to the position dependence of the 
diffuse X-ray emission in the galaxy, the PSF extent, the 
exposure, and CCD efficiency. Therefore it is necessary 
to carry out the incompleteness correction for the XLF. 
The D 25 region of NGC 3115 has a semi-major axis of 
a = 3.62' (10.2 kpc), a semi-minor axis of b = 1.23' (3.5 
kpc ) and a position angle of 40° (Ide Vaucouleurs et al.l 
ll99lh . To limit the incompleteness effects and the cosmic 
X-ray background (CXB) contribution, we defined our 
study field for the XLF of field LMXBs as the region 
inside D 25 . Further considering that the central region 
is very crowded and has strong diffuse X-ray emission, 
we excluded the central a = 10 " elliptical region (the 
eccentricity and position angle follow the D 25 ellipse) for 
all XLFs throughout the paper. 

We calculated the incompleteness function K(L), the 
fraction of pixels weighted by the assumed spatial distri¬ 
bution of sources, in which a source with the luminosity L 
or higher wou ld be detected, using th e backward correc¬ 
tion method dKim fe Fabbianol l2003fl . In this method, 
sources are simulated with MARX and added one by 
one to the real observed image, which is then checked 
to see whether each one could be detected with uiavde- 
tect. The source spectral shape was assumed to be a PL 
with a photon index of 1.7 and the Galactic absorption. 
The simulations were carried out for a series of lumi¬ 
nosities with an increasing factor of 1.1 and 1.21 below 
and above the 90% completeness luminosity, respectively. 
The positions of simulated sources for each luminosity 
were specified as follows. The D 25 ellipse was divided 
into elliptical annuli with a series of ellipses that have 
the eccentricity and the position angle following D 25 and 
semi-major axis a for the ith (i = 0 to 56) ellipse as¬ 
suming ai = ai -1 + Sa * 1.04® -1 , where a 0 = 10" and 
Sa = 1". Each elliptical annulus was then divided into 
80 cells with an equal area, and the simulated source po¬ 
sition was specified at the center of each cell, with 4560 
in total. In some cases the simulated sources coincided 
with the real sources, and we assumed that the simulated 
sources were detected by wavdetect only if the simulated 
sources dominate the flux over the real sources, which 
is to take into account the source confusion effect. We 
calculated K(L) for the CXB sources and field LMXBs 
separately, because the CXB sources have a flat distri¬ 
bution and the field LMXBs are expected to follow the 
Ag-band light dGilfanovI I2004H . f or which we used the 
2MASS Large Galaxy Atlas data (jJarrett et al.ll2003h . 

We also calculated K(L) for GC LMXBs. We used all 
the GCs detected in the optical (Section 12.311 as the par¬ 
ent spatial distribution of GC LMXBs and assumed that 
they have equal probability of hosting an LMXB. To limit 
the incompleteness effects and spurious rate of the GC 
LMXB identification, our study field for the XLF of GC 


LMXBs is set to be the HST /ACS field of view (FOV). 
The HST /ACS FOV reached ~1.3 £> 25 and ~2.5 D 25 in 
the major-axis and minor-axis directions of the D 25 el¬ 
lipse, respectively (Figure CD- 
The differential XLF of LMXBs in a given regi on can 
be calculated as follows (refer to, e.g. JVoss et al.l[200^ 1: 


^-Wmxb _ 1 f dN obs 

dL A'lmxb(A) \ dL 


- 4itD 2 K C xb{L) 


dNcxB \ 
dL ) 
( 1 ) 


where N obs is the total number of observed sources 
and D is the distance to NGC 3115. The quantity 
4nD 2 dNcxB/dL is dNcxB/dS, the log (N) log(S') distri¬ 
bution of the CXB sources. We used the full band (0.5- 
10 keV) log (TV)-log (5*) d istribution of CXB sources from 
iGeorgakakis et ahl (120081) . with their 0.5-10 keV flux con¬ 
verted to our 0.5-7 keV band assuming a PL spectrum 
with a photon index of 1.4. For the XLF of field LMXBs 
only, we filtered out GC LMXBs and had the CXB con¬ 
tribution estimated as described above. For the XLF of 
GC LMXBs, the CXB contribution was not corrected 
because it is negligibly small. 

We did n ot correct XLFs f or the HMXB contribution. 
Following IMineo et all (I2012H , we estimated the star for¬ 
mation rate in NGC 3115 to be 0.07 M@ yr -1 (see their 
Equation 9, which is based on the UV and IR emission). 
Based on their XLF for HMXBs (their Equation 18), we 
can estimate the number of HMXBs in NGC 3115 above 
Li; m to be 2.3, which is one order of magnitude less than 
the CXB contribution and is thus negligibly small. 

To compare with previous studies, some XLFs pre¬ 
sented in this study will be divided by (thus normalized 
to) the stellar mass enclosed in our study region of the 
field LMXB XLF (i.e., within £>2 5 and outsi d e the cen¬ 
tral a = 10" ellipse). Following IZhang et al.1 (l2012h . we 
used the A' s -band luminosity and estimated the stellar 
mass in our study region to be 6.31 x 10 10 M 0 (the total 
stellar mass within D 25 is 7.83 x 10 10 M 0 ). 


2.3. Multiwavelength cross-correlation 

We cross-correlated our X-ray sources with op¬ 
tical sources from HST /ACS mosaic imaging and 
Subaru/ Suprime-Cam imaging to search for the GC 
LMXBs. I Jennings et al.l (1201411 compiled 360 GC candi¬ 
dates from HST /ACS mosaic imaging and an a dditional 
421 from Subaru/Suprime-Cam imaging (lArnold et all 
I20T1 . Before cross-correlation, we first carried out 
absolute astrometry correction on X-ray sources by 
cross-correlating their positions with the 360 GC candi¬ 
dates from HST /ACS mosaic imaging, who se astrometry 
was r egistered to the USNO-Bl.O Catalog dMonet et al.l 
I2003T) . We only used X-ray sources detected at > 6 a 
significance and with off-axis angles < 6 ' (the limit of 
HST /ACS FOV) in the cross-correlation. We found 30 
matches with a median separation residual of 0.06". 

We searched for the HST /ACS GC counterparts to our 
X-ray sources using the 99.73% (i.e., 3cr) positional un¬ 
certainty that combines both X-ray and optical compo¬ 
nents. For the HST /ACS sources, we assumed the half 
light radius as the 1-a positional uncertainty. We also 
included a systematic uncertainty which was assumed to 
be 0.05" (lcr, in both R.A. and Deck) based on the above 
matches in the absolute astrometry correction. This sys¬ 
tematic uncertainty is probably overestimated, but it is 
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so small that the number of GC matches remained the 
same even if we did not include this systematic uncer¬ 
tainty. The offsets of all matches (37 in total) are < 1" 
(only 3 have offsets > 0.2") and have a median of 0.07". 
To estimate the spurious rate, we rotated the HST /ACS 
field by ±10°, 180°±10°, and 180° around the galaxy 
center and carried out the cross-correlation in the same 
way and found the spurious rate to be about 3%. _ 

The HST /ACS GC distribution from iJennines et al.l 
l20Tl decreases sharply within ~O. 25 .D 25 (see Sec¬ 
tion liTID . Therefore their GC list is probably fairly in¬ 
complete in this region owing to strong stellar light. We 
tried to match ou r X-ray sources with the sources that 
were detected by I Jennings et al.l (12014T ) but not classi¬ 
fied as GCs, and we found three extra matches (S12, S53 
and S79 in Paper I) within O. 25 D 25 . One more source 
(S65) seems to have an optical ma tch from our visual in- 
spection but it is not detected by Uennings et, al.l (|2014T ) 
due to its being too close to the galaxy center (5.5", 
within the a = 10" elliptical exclusion region). We clas¬ 
sify these four sources as GC LMXB candidates. They 
are all bright (>10 37 erg s -1 ) and are not expected to 
be CXB sources (the expected CXB source number at 
this luminosity is 0.5 and is much smaller if only CXB 
sources with bright optical counterparts are considered). 
Another source (S92 in Paper I) outside 0.25-D 2 5 also 
has an optic al match , but it was classified as a star by 
LTennings et al.l ( 2014 1 due to the measurement of a ra¬ 
dial velocity (238 km s -1 ) much lower than the threshold 
of 350 km s -1 that they adopted to define GCs. Con¬ 
sidering that the size and color of this optical match are 
consistent with typical GCs, we treat it as a GC LMXB 
candidate too. We did not include the above five LMXBs 
in either field or GC XLFs (although S65 is outside the 
study region and would not be included anyway). 

In our search for the GC optical counterparts to our 
sources detected only in the Subaru/Suprime-Cam imag¬ 
ing (i.e., not in the HST /ACS imaging), we also used the 
99.73% positional uncertainty. The 1-cr positional uncer¬ 
tainty of the optical sources was assumed to be 0.1" in 
both R.A. and Deck The 1-cr systematic uncertainty was 
assumed to be 0.1" in both R.A. and Deck (Paper I). To 
limit the spurious rate, we have a maximum searching 
radius of 2". The spurious rate was estimated to be 5% 
(Paper I). Only eight Subaru/Suprime-Cam GC matches 
were found. They are used only for the study of the spa¬ 
tial distribution of GC LMXBs in Section l3Jl but not for 
the study of XLFs. 

3. RESULTS 

3.1. Spatial Distribution 

Figure Q] shows the Chandra X-ray image of NGC 
3115. From the merged observation, we detected 145 
X-ray sources above the 50% completeness luminosity 
L\\ m = 1.21 x 10 36 erg s _1 (Section [3.2D within D 25 and 
outside the central a = 10" ellipse (23 within this exclu¬ 
sion region). The expected CXB source number is 26.5. 

Figure [2] shows the cumulative radial distributions of 
different classes of objects detected from the merged 
observation. Due to the high inclination (i = 86°, 
iCanaccioli et alJ 1 19871) of NGC 3115, we plot the dis¬ 
tributions with respect to a/R 25 , where a is the angular 
separation between the source and the galaxy center and 



a /^25 

Fig. 2. — The radial distribution of observed sources excluding 
GC LMXBs (red solid line), compared to the model (thick blue 
dotted line) composed of field LMXBs (FLMXB) and CXB sources 
(thin dotted line). Also plotted are the radial distributions of GC 
LMXBs (red dot-dashed line) and the optical GCs (dashed line, 
divided by a factor of 10). The cumulation starts at a/i ?25 = 
0.046 (i.e., excluding the central a = 10 // elliptical region) and 
ends at 2 D 25 . To reduce the incompleteness effects, only sources 
with Lx > 4.0 x 10 36 erg s -1 are used. 



Fig. 3. — Cumulative XLFs from the merged observation (black 
solid line) and single observations 2040 (red dotted line), 12095 
(green dashed line), and 13820 (blue dot-dashed line). The grey 
solid line is also from the merged observation but using sources 
with Vvar <5.0. The grey dot-dot-dashed line uses the maximum 
luminosity of each source. The distributions are not corrected for 
incompleteness or the CXB contribution. 






























5 



Fig. 4.— The incompleteness functions for field LMXBs (red 
solid line), GC LMXBs (green dashed line) and CXB sources (blue 
dotted line). 

i ?25 is the elliptical radius of the £>25 isophotal ellipse in 
the direction from the galaxy center to the source. The 
cumulation starts at a/R 2 5 = 0.046 because we have 
excluded the central a = 10" elliptical region in the cal¬ 
culation of XLFs. To reduce the incompleteness effects, 
we used only sources above 4 x 10 36 erg s _1 , which is 
the 82% completeness luminosity for CXB sources and 
the 94% completeness luminosity for field LMXBs within 
2£>25 (but outside the central a = 10" elliptical region). 
The non-GC X-ray sources, expected to consist of field 
LMXBs and CXB sources mostly (red solid line), and 
GC LMXBs (red dot-dashed line) are plotted separately. 
The thick blue dotted line models the spatial distribution 
of non-GC X-ray sources using two components, one for 
CXB sources and the other for field LMXBs, which was 
assumed to follow the IR light in the I\ s band. The nor¬ 
malization of the field LMXB component was determined 
so as to give the same number of sources, after adding 
the CXB component, as observed within 2D 2 $. The ob¬ 
served distribution roughly follows this model, with the 
Kolmogorov-Snrirnov (K-S) test giving a probability of 
30%. The radial distribution of GC LMXBs seems to 
approximately follow that of GCs detected in the optical 
as well, with the K-S test giving a probability of 62%. 

Figure [2] shows a dramatic difference between the spa¬ 
tial distributions of field and GC LMXBs. The field 
LMXBs tend to cluster toward the galaxy center, while 
GC LMXBs tend to be more spread out. Only 7 out of 
360 AST/ACS GCs are within 0.2Z?25- However, we note 
that the GC detection is most probably fairly incomplete 
near the galaxy center due to strong stellar light. Some 
GCs could also be destructed near the galaxy center due 
to mass segregation. 

There are 50 non-GC X-ray sources above 4 x 10 36 erg 
s _1 observed between D 2 5 and 2 D 2 5 , while the expected 
number of CXB sources is 40.2 after incompleteness cor¬ 
rection and the expected number of field LMXBs are 1.5 
based on the IR light in the K s band. The above 50 
sources include two special sources: S109 and S179 (Pa¬ 
per I). The former is a supersoft X-ray source (SSS) at 


table 2 

The maximum likelihood fits to XLFs of different populations 

USING A BROKEN POWERLAW 


Population 

iVtot 

NcxB 

ai 

Oi 2 

L b 

K 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Total 

145 

26.5 

1 03+ 010 

l.UO — o.14 

i -) 9+0.3 
Z - Z -0.4 

55l^ 

30+ 

1 1 q + 0.08 

i ' 1 °Q,Q8 

> 6.7 

164±| 3 

36+ 

Field 

114 

26.5 

1 io" l_u - ±z 

11U -0.13 

o 4~ l u '° 

z -*-0.5 

57+ 

23+ 

-1 9-1 +0.09 
± -^ ± - 0.09 

> 7.0 

198+13 

32+ 

Fieldi n 

44 

1.0 

0.84±°; 2 3= 

0.6 

33+5 

9+ 

Fieldout 

70 

25.5 

i-20 ++ 

> 9.6 

1661" 

15+ 

GCacs 

36 

0 

0 lO" 1 " 0 ^ 
u - ±u -0.62 

1 6" t uz 

J -- u -0.2 

H + 

1-8+.2 

GCd25 

27 

0 

0 41+0-39 

u -^ x -0.56 

1 «+ 0 - 2 

12+ 4 

9 c+2-8 

Note. — 

Columns 

are as 

follows. (1) 

i The population, 

(2) the to- 


tal number of sources observed above Ln m — 1.21 X 10 36 erg s —1 , (3) 
the expected obse r ved n umber of CXB sources above based on 

IGeorgakakis et al.| d2008l) , (4) initial slope, (5) second slope, (6) break 
luminosity in units of 10 36 erg s —1 , (7) the normalization (not nor¬ 
malized by the stellar mass). All errors and lower limits are at the 
lcr level. Population descriptors: “total”: all LMXBs within (0.046— 
1.0)^25, “field”: all LMXBs within (0.046—1.0)Z)25 but excluding all 
27 GC LMXBs and four GC LMXB candidates in the region, “fieldin’ 5 : 
similar to “field” but only within (0.046—0.2)L>25> “field out ”: similar 
to “field” but only within (0.2—1.0)^25 5 GCacs ; all HST /ACS GC 
LMXBs, GCd 25 : all HST /ACS GC LMXBs within D 25 - For the “to¬ 
tal” and “field” XLFs, we give both the low-break and high-break 
solutions. 

I. 53 D 25 . The latter is a transient with 0.5-7 keV long¬ 
term variability factor War > 20 and relatively soft X-ray 
spectra (classified as a BH X-ray binary candidate in Pa¬ 
per I), and it is at 1.04D 2 5. Therefore these two sources 
are most probably in NGC 3115, instead of being CXB 
sources. Then we have 48 left, which is 19% more than 
the expected number of CXB sources. If we concentrate 
on the region between 1.3Z?25 and 2D 2 $ (there is little IR 
light outside 1.3£>25 from the galaxy), we have 33 non- 
GC X-ray sources above 4 x 10 36 erg s _1 , excluding S109, 
and this number is very close to the expected number of 
CXB sources (30.6, incompleteness corrected). 

3.2. The XLF of All LMXBs 

Because a large fraction of our sources are variable, 
we first check whether there is any difference in XLFs 
between observations. We plot the observed cumulative 
XLFs for the merged observation and the longest ob¬ 
servation in each of the three epochs in Figure [3] We 
compared their XLFs using the K-S test. We focused on 
luminosities above the 90% completeness limit, which are 
3.3 x 10 37 erg s _1 , 1.6 x 10 37 erg s _1 , and 8.6x 10 36 erg s _1 
for observations 2040, 12095 and 13820, respectively. We 
found that the XLFs obtained in these individual obser¬ 
vations are consistent with the same distribution as that 
obtained in the merged observation with the K-S prob- 
ability of 48%, 98%, and 82%, respectively. lYoss et all 
(|2009D obtained a similar conclusion for Centaurus A. 
We also constructed the XLF using less variable sources 
(the 0.5-7 keV flux long-term variability War < 5, Paper 
1) in the merged observation (gray solid line in Figure [3]) 
and compared it with the XLF using all sources. Above 
the luminosity at the 90% completeness level (2.7 x 10 36 
erg s _1 ), the K-S test also indicates no obvious difference 
between them, with the probability of 100%. Finally, we 
constructed the XLF using the maximum luminosity of 
each source (gray dot-dot-dashed line in Figure[3]). When 
we compared it with the XLF from the merged observa- 
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Fig. 5. — Incompleteness corrected and CXB contribution subtracted XLFs for all LMXBs (left panels), field LMXBs (middle panels) 
and GC LMXBs (right panels). We plot both the cumulative (upper panels, the red solid line with the shaded area representing la 
Poissonian uncertainty) and differential forms (lower panels). Our fits with a broken PL are shown as green solid lines (the low-break 
solution) and purple dot-dashed lines (the high-break solution, when present, see text). We note that the high-break solutions have a large 
uncertainty in the second index, whose lower limit is given in Table [2] The dashed blue lines in the left panels are the fit to the average 
XLF of 20 early-type galaxies bv I Zhang et al.l H2012T) . with the normalization decreased by 24%. The black dotted line is the expected CXB 
distribution. 


tion and limited to L > 2.7 x 10 36 erg s _1 , the K-S test 
gave a probability of 37%. From all the above compari¬ 
son, we see no significant effect of the source long-term 
variability on the XLFs. Therefore, we will use the mean 
luminosities of the sources detected from the merged ob¬ 
servations for the XLFs hereafter. 

The incompleteness functions for (field) LMXBs and 
CXB sources are shown in Figure 0] While they ap¬ 
pear to be similar to each other, this is a coincidence, 
as they can be very different if different regions other 
than the whole D 25 (excluding the central a = 10" el¬ 


lipse) are used. The luminosity corresponding to the 50% 
completeness level is about L\ im = 1.21 x 10 36 erg s _1 , 
which is much lower than typical values of ~ 10 37 erg s -1 
achieved for other galaxies by Chandra (e.g. dZhang et ahl 
I2012T) . Above this luminosity limit, the sources detected 
in the merged observation within D 25 are all >2.9 ct. 

The incompleteness corrected and CXB contribution 
subtracted XLF above Ln m for all LMXBs is plotted in 
the left panels in Figure [5] The XLF is steep down to 
a break around 5 x 10 37 erg s _1 , below which the XLF 
flattens clearly. We fitted the differential form of XLF 
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Fig. 6.— Incompleteness corrected and CXB contribution sub¬ 
tracted XLFs of field LMXBs within a/ifes = 0.046-0.2 (filled cir¬ 
cles) and those within a/i ?25 = 0.2-1.0 (red open circles). Their 
best-fitting broken PL is shown as a dotted line and a red solid 
line, respectively. 



Fig. 7. — The color-magnitude diagram of HST/ ACS GCs, with 
those containing LMXBs enclosed with squares (Lx within (0.8— 
1.3)xlO 37 erg s —1 ) or circles (others), whose size is proportional to 
the logarithm of the LM XB l umi nosity. The re d dotted line at g — 
z = 1.13 mag was used bv l Jennings et all <120141) as the dividing line 
between blue/metal-poor and the red/metal-rich subpopulations. 

with a small Lx bin size of S log(Lx) = 0.02 using the C 
statistic (it is maximum likelihood-based) in the Xspec 
fitting package (lArnaudlll996l) . In this way, the observed 
XLF was used in the fit, but the fitting model was mod¬ 
ified by the incompleteness function (through a response 
file). The model that we adopted is a broken PL: 

dN _ \KL^\ L 36 < L b ( ) 

dL 36 {KL^-^L^ 2 , L 3 6 > L b ’ U 

where L 36 = Lx/(10 36 ergs -1 ). We obtained a\ = 


1.03tg;^, a 2 = 2.2±g;|, and L b = 55t^ (Tableland 
F igure [5] (green solid line)). 

IZhang et all (|2012f) fitted the average XLF of 20 early- 
type g alaxies with the template introduced by iGilfanovl 
(2004|), which is essentially a double broken PL. Be¬ 
cause we do not have enough statistics above the sec¬ 
ond break that they obtained (6 x 10 38 erg s -1 ), we 
do not need to introduce the second break to fit our 
XLF. Our fit is fully consistent with that obtained by 
IZhang et all (|2012f l. who reported a\ = 1.02 ± 0.08, 
a 2 = 2.06 ± 0.06, and L b = 54.6 ± 4.0. The main dif¬ 
ference is the normalization A', with ours being about 
76% of that of IZhang et alJ (|2012f l. after being normal¬ 
ized by stellar mass (our I\ = 4.6 ± 1.3 per 10 10 Mg 
versus their K = 6.0 ± 1.7 per I Q 10 M^) . This is c onsis- 
tent with the result obtained bv IZliang et all (12 0121 ) that 
NGC 3115 has a relatively low number of LMXBs per 
unit stellar mass compared with other galaxies studied 
by them (NGC 3115 was included in their galaxy sample 
with only the first three observations analyzed). Their fit 
with the normalization decreased by 24% is also shown 
in Figure [5] (blue dashed line). 

We note that we also found a high-break fit (Table 0 
purple dot-dashed line in Figure 0 that gives a higher 
break luminosity ((1.6 ± 0.1) x 10 38 erg s -1 ) and a much 
steeper second slope (> 6.4, 1 -a lower limit) and has a 
C statistic only larger than the above low-break fit by 
1.1. The initial slope of the high-break fit (1.13 ± 0.08) 
is slightly higher than that of the low-break fit. 

3.3. The XLF of LMXBs in the Field 

The incompleteness corrected and CXB contribution 
subtracted XLF above L\[ m for LMXBs in the held is 
plotted in the middle panels in Figure 0 The results 
of our fit with a broken PL are given in Table [2] As for 
the XLF of all LMXBs, we also find a low-break fit and a 
high-break fit to the XLF of held LMXBs. The low-break 
ht has a C statistic higher than the high-break ht by only 
0.3. The parameters of both hts are very similar to the 
corresponding hts to the XLF of all LMXBs (Section l3.2L 
probably due to relatively few GC LMXBs in the total 
sample. Figure 0 shows both the low-break ht (green 
solid line) and the high-break ht (purple dot-dashed line) 
to the held LMXB XLF. 

To check whether there is spatial variation of the XLF 
of held LMXBs, we divided our study held into an inner 
((0.046-0. 2)D 25 ) and an outer ((0.2-1.0)D 2 .5) region (the 
boundary was chosen to ensure enough statistics in both 
regions) and created two corresponding XLFs. They are 
shown in Figure 0 It is nontrivial to use the K-S test 
to compare these two XLFs due to their different CXB 
and incompleteness corrections. Therefore, we also fit¬ 
ted them with a broken PL for comparison. The htting 
results are shown in the hgure and given in Table 0 The 
XLF of held LMXBs in the inner region seems to be hat¬ 
ter at low luminosities, and the break seems to be at a 
lower luminosity than that in the outer region. We see no 
clear degeneracy in the ht any more, but the ht to the 
XLF from the inner region is similar to the low-break 
ht to the XLF from the whole region, and the ht to the 
XLF from the outer region is similar to the high-break ht 
to the XLF from the whole region, especially the break 
and second slope. The best-fitting an, L b , and 0:2 of the 
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XLFs in the inner and outer regions differ at the 1.5er, 
1.8a, and 1.5a confidence levels, respectively. Such dif¬ 
ferences are marginally significant, and we will discuss 
the possible origin in Section 0 

3.4. The XLF of LMXBs in GCs 

The XLF for the LMXBs detected in HST/ ACS GCs 
are shown in the right panels in Figure [5] It is in¬ 
completeness corrected using the incompleteness func¬ 
tion (green dashed line) shown in Figure [I] The XLF of 
GC LMXBs seems flatter than that of field LMXBs. The 
fit with a broken PL is given in Table [2] We obtained 
on = 0.10±Q g2, «2 = 1-6 ± 0.2, and Lb = lll^, which 
are different from those from the field LMXBs (the low- 
break fit) by 2.5cr, 1.6a and 2.3d, respectively. Our re¬ 
sults agree with previous finding that the XLF of th e field 
LMXB is steeper than that of GC LM XBs dVoss et al.l 
120091 [Kim et al.l[2009HZhang et alJ[2011l h though the dif¬ 
ference is not very significant in our case, owing to the 
relatively few GCs in our sample. 

Previous studies of the XLF of GC LMXBs often just 
used sources within the D 25 ellipse. We also obtained 
such an XLF for GC LMXBs, and the fitting results with 
a broken PL are given in Table [2] We obtained (\\, 02, 
and Lb different from those from the field LMXBs (the 
low-break fit) by 1.5a, 1.5a, and 1.1a, respectively. Due 
to fewer sources used, these differences are less significant 
than those reported above using all HST/ ACS GCs. 

3.5. Optical Properties of GC LMXBs 

Figure [7] plots the color-magnitude diagram of 
HST/ ACS GCs. To indicate the presence of LMXBs and 
their luminosity we circle the 37 GCs containing LMXBs, 
with the size of the circle proportional to the logarithm of 
the LMXB luminosity, but for 13 GCs whose Lx clusters 
within (0.8-1.4)x 10 3 ' erg s -1 , they were enclosed with 
squares instead. We note that all our GC LMXBs have 
Lx > 1.13 x 10 36 erg s _1 . The median ( g — z) colors of 
GCs with LMXBs and without LMXBs are 1.23 and 1.06, 
respectively. Based on the nonparametric Wilcoxon rank 
sum test, the difference is at a significance level of 2.7a. 
Concentrating on luminous LMXBs with Lx > 10 37 erg 
s _1 (24 in total), we found that the median (g — z ) col¬ 
ors of GCs with luminous LMXBs and without luminous 
LMXBs are 1.31 and 1.05, respectively, corresponding 
to a Wilcoxon rank sum difference of 3.5a. If we fol¬ 
low |jennmgijjtjd] (12014T) and use g — z = 1.13 mag as 
the dividing line (dotted line in Figure 0 between the 
red/metal-rich and blue/metal-poor subpopulations, we 
find that 23 out of 169 (i.e., 13.6%) red GCs contain 
LMXBs, while there are 14 out of 191 (i.e., 7.3%) blue 
GCs containing LMXBs. Thus, the fraction of red GCs 
hosting LMXBs is about twice of that of blue GCs host¬ 
ing LMXBs. Concentrating on luminous LMXBs with 
L x > 10 37 erg s” 1 , we find 18 (i.e., 10.7%) red GCs 
and 6 (i.e., 3.1%) blue GCs hosting luminous LMXBs. 
The former fraction is 3.5 times of the latter, which is 
consistent with previous studies using a sample of galax¬ 
ies with limiting X - ray luminosity arou nd 10 37 erg s _1 
(|Kundu et al.ll2007t iSivakoff et alJl2007h . 

Figure 0 also shows that GCs hosting LMXBs tend to 
be bright/massive, as found previously fo r many galaxies 
(|Sivakoff et al.ll2007t iKundu et al.l I2007L including NGC 


3115, but using only the Chandra observation in 2001). 
The median M z is —9.62 for GCs with LMXBs and 
is —8.23 for GCs without LMXBs, corresponding to a 
Wilcoxon rank sum difference of 6.1a. Similar results 
can be obtained if we just focus on luminous LMXBs 
with Lx > 10 37 erg s~ 3 , with the median M z of —9.48 
for GCs with luminous LMXBs and —8.29 for GCs with¬ 
out luminous LMXBs (the Wilcoxon rank sum difference 
is 4.4a). Separating the GC subpopulations and concen¬ 
trating on bright GCs with M z < —9.0 mag, we find 16 
out of 53 (i.e., 30%) red GCs and 10 out of 36 (i.e., 28%) 
blue GCs hosting LMXBs. For even brighter GCs with 
M z < —10, we find 10 out of 14 (i.e., 71%) red GCs 
and 3 out of 8 (i.e., 38%) blue GCs containing LMXBs. 
These fractions are significantly higher than those ob¬ 
tained above for all GCs (i.e., 13.6% and 7.3% for red 
and blue GCs, respectively). 

4. DISCUSSION 

4.1. The Correction of CXB Contribution in XLFs 

We have obtained the XLFs of LMXBs in NGC 3115 
down to Li; m rs 10 36 erg s _1 , which has only been 
achieved fo r one other old galaxy, i.e., Centaurus A 
(jVoss et al.l I2009T ). We have carried out careful correc¬ 
tions to the XLFs to account for the incompleteness ef¬ 
fects and CXB contribution. We found no large discrep¬ 
ancy bet ween the CXB density in o ur field and that esti¬ 
mated bv lGeorgakakis et all (12008t h who used data from 
six large Chandra surveys. Even assuming a possible 
20% enhancement of the CXB density in our field (Sec¬ 
tion [30, we found no noticeable effect on the XLFs. In 
some studies, the CXB contribution was taken into ac¬ 
count by directly excluding CXB so urces id e ntified from 
the optical cross-correlation fe.g.. lKim et al.l[2009ih We 
did not show the results using this method because our 
HST imaging is not deep enough. However, we also 
tested this method by excluding the AGNs that we iden¬ 
tified in Paper I (about 50% of the expected number) 
and obtained XLFs very similar to the ones that we have 
shown. This is mainly because of the high inclination of 
NGC 3115 so that the CXB contribution is less signifi¬ 
cant for this galaxy than other typical ones. 

4.2. Caveats on the XLF of Field LMXBs and Physical 
Implications 

Our XLF of field LMXBs can be fitted with parameters 
typically seen in the literature for other old normal galax¬ 
ies or the bulge of spiral galaxies, with a possible break 
around 5.7 x 10 37 erg s _1 . Such a break has been at¬ 
tributed to different mechanisms of removal of orbital an¬ 
gular momentum (mag netized ste llar wind versu s gravi¬ 
tational wave emission, iPostnov &; Kuranovll2005l ) or dif¬ 
ferent ty£es_of_donor_stars_(giants versus main-sequence 
stars. iRevnivtsev et al.ll2011f) in the high and low lumi¬ 
nosities. 

However, the interpretation of our XLF of field LMXBs 
is complicated by the presence of a degenerate solution 
with a high break at about 2 x 10 38 erg s _1 . This degen¬ 
eracy appears to be associated with a spatial variation 
of the XLF, with the XLF in the inner region ((0.046- 
0.2)1)25) being flatter and having a lower break (3 x 10 3 ' 
erg s _1 ) than that in the outer region ((0.2-1.0)D 2 5), 
which has a break around 1.7 x 10 38 erg s _1 . Therefore, 
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the XLF of field LMXBs in the inner region is closer to 
the XLF of GC LMXBs. One possible cause for this is 
that our field LMXB sample in the inner region could 
include some GC LMXBs that we cannot identify due 
to significant incompleteness effects near the galaxy cen¬ 
ter that limit our ability to detect GCs in the optical. 
In the outer region we found 25 LMXBs from 219 GCs 
and 44.5 in the field (after excluding the CXB contribu¬ 
tion) above Ln m , but only 2 GC LMXBs from 7 GCs 
and 43.0 field ones in the inner region. To have the 
same ratio of GC LMXBs to field LMXBs, we would have 
missed 22 GC LMXBs in the inner region. However, it is 
well known that the spatial distribution of GCs is more 
extended than the stellar light, as can also be seen in 
Figure [2] for the outer region. The distribution of GCs 
within (0.3-1.3)1)25, if fitted with a PL, is found to follow 
- \/a/R 25 , which would indicate that in the inner region 
there should be 93 GCs. Assuming the same detection 
rate of LMXBs (25/219 = 11.9%) as in the outer region, 
we would have missed 6.6 GC LMXBs, given that we 
have detected 2 GC LMXBs and 2 candidates. Thus, 
the number of GC LMXBs that we missed in the inner 
region is probably small, compared with the total number 
of sources observed (43.0 after subtracting the CXB con¬ 
tribution), and their effect on the XLF of field LMXBs 
in the inner region should be insignificant. 

Alternatively, the spatial variation of the XLF of field 
LMXBs might be real and can be explained if the field 
LMXBs in the inner region have a dynamical origin sim¬ 
ilar to GC LMXBs. There are two scenarios: one is 
the dynamical formation of LMXBs in the dense stel¬ 
lar environment near galaxy nuclei, and the other is the 
destruction of GCs that drift toward the galaxy center 
due to mass segregation, leaving behind the remnant 
LMXBs. The former was argued to be the dominant 
mechanism to account for the high specific frequency of 
X-ray sources, per unit stellar mass (following the pi de¬ 
penden ce on the stellar d e nsity), near the center (<F) of 
M31 bv lVoss fc Gilfanovl (l2007aM . However, we do not 
see increasing high specific frequency of (non-GC) X-ray 
sources at the very center, compared with the outer re¬ 
gion (Figure [2]). The stellar de nsity of M31 is around 3 0 
M 0 pc -3 at V from the center (IVoss &; Gilfanovll2007blb 
Based on the stellar density model bv lEmsellem et all 
Gael), we expect that NGC 3115 reaches a similar stellar 
density at 8" in the major-axis direction. Therefore 
the former mechanism is probably still not significant in 
our inner region, which excludes the central a = 10 ,/ ellip¬ 
tical region. There is large uncertainty in estimating the 
level of the second mechanism. The specific frequency of 
(non-GC) X-ray sources seems to peak around a/i?25 ~ 
0.1-0.2. If some part of it is due to the second mecha¬ 
nism, the remnant LMXBs should gain some momentum 
to reach this region during the destruction of GCs, or the 
destruction should be able to occur there. 

Considering the possible contamination of the LMXBs 
from undetected or disrupted GCs in the inner region, the 
XLF in the outer region is probably a better represen¬ 
tation of the characteristics of primordial field LMXBs. 
Its break at about 1.7 x 10 38 erg s _1 (Table [2]) seems 
somewhat higher (at the 90% confidence level) than 
typical values of around 5 x 10 37 erg s" 1 reported in 
other studies that normally used most of the field (that 


is, no differentiation between the inner and outer re¬ 
gions! dGilfanovI 12004 1 Voss et al.ll2009HKim et all [20091 
IZhang et all l2012f) . iZhang et al.l T201.3l) obtained XLFs 

combining 20 early-type galaxies for the inner and outer 
regions separately. Their inner and outer regions were 
defined as (0.2-3)r e and (4-10)7' e , respectively, where r e 
is the Ag-band half-light radius. Our inner and outer 
regions for NGC 3115 are approximately (0.3-1.3)r e and 
(1.3-6.5)r e , respectively. Visually it appears that the 
break of the XLF in their outer region is higher than 
the XLF in their inner region (see their Figure 4). How¬ 
ever, they did not carry out the fit, and the significance 
of this variation is not clear. Moreover, they did not 
exclude GC LMXBs, making it difficult to compare di¬ 
rectly. In the future, more galaxies should be used to 
investigate the XLF of field LMXBs from the outer re¬ 
gion to check whether the high break that we observed 
in the XLF of field LMXBs in the outer region of NGC 
3115 is universal or due to statistical fluctuation. If it 
is real, the best explanation is probably the Eddington 
limit of NS LMXBs. In Paper I, we have shown that 
most of our bright LMXBs (above several 10 36 erg s _1 ) 
are NS LMXBs in the soft state. 

In summary, the flatter XLF of field LMXBs in the 
inner region compared to that in the outer region is un¬ 
likely due to dynamically formed LMXBs in the dense 
stellar environment near the galaxy nucleus, but could be 
due to contamination of LMXBs from undetected and/or 
(more likely) disrupted GCs in the inner region. The 
field LMXBs in the outer region are more likely to be 
primordial. The break of their XLF could be due to the 
Eddington limit of NS LMXBs, agreeing with our find¬ 
ing in Paper I that most of our bright sources are NS 
LMXBs in the soft state. 


4.3. Caveats on the XLF of GC LMXBs and Physical 
Implications 

Considering the large detection rate of LMXBs (~71%, 
Section [33]) in the most metal-rich and the most massive 
GCs, some of these GCs probably in fact host multiple 
LMXBs which cann ot be resolved by Chandra. We fol¬ 
low the method of iSivakoff et all (2 0071 ) to study such 
source blending effects. They found the dependence of 
the expected number At (assuming Poisson statistics) of 
LMXBs per GC on the GC properties to be: 


At — A 


M 

10 6 M q 


1.237 

l(]0.9(g-z) 



( 3 ) 


where the GC mass is M = 1.45 x 1 Mq 

= 4.512) and the half-mass radius is rh,M = 
rh x 10°- 17 [G -Z ) -1 - 22 ]. To match our observation of 37 
GCs hosting LMXBs, the normalization, A, should be at 
least 0.16, without taking into account the incomplete¬ 
ness effects. The expected median M z is then —9.36, 
sligh tly f ainter than the observed value (i.e. —9.62, Sec¬ 
tion [33]). Concentrating on the region of M z < —10.0 
and (g — z) > 1.13 where the LMXB detection rate is 
the highest, with 10 out of 14 GCs observed to host 
LMXBs (Section [331) . the expected number of GCs host¬ 
ing LMXBs from Equation [3] with A = 0.16 is 8.0, with 
the number of GCs expected to host multiple LMXBs 
A/nuiti = 3.5. The total expected number of LMXBs 
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is Almxb = 13.5. To match the observed value of 10 
GCs hosting LMXBs, A should be 0.25, resulting in 
A r mu iti = 5.8 and TVlmxb = 21.6 in this region. A similar 
exercise for the rest of the region (i.e., M~ > —10.0 or 
(g — z) < 1.13), we find that A = 0.14 could match the 
observed 27 GCs hosting LMXBs, with iV mu i t i = 3.4 and 
Almxb = 31.1. 

GCs hosting multiple LMXBs are expected to be more 
lumino us a nd show less long-term variability in X-rays 
dKundu et al.1 12007t h In the dense LMXB region of 
M z < —10.0 and (g — z) > 1.13, there are four with 
the maximum 0.5-7 keV luminosity > 5 x 10 3 ' erg s -1 . 
They have 0.5-7 keV long-term luminosity variability of 
V = 1.4, 1.4, 2.1, and 14.9. Among the field LMXBs 
within (0.046-1.0)I?25, there are 27 with maximum 0.5-7 
keV luminosity > 5 x 10 37 erg s _1 . Their median vari¬ 
ability is 2.3, which is not significantly larger than that 
found above for the GC LMXBs with high likelihood of 
blending. Thus we cannot confirm any effect of blending 
on the variability. 

To investigate the source blending effects on the XLF, 
we carried out Monte Carlo simulations using Equation[3] 
with A = 0.25 if M z < —10.0 and (g — z) > 1.13 and 
A = 0.14 elsewhere, as obtained above. When a GC was 
simulated to host multiple LMXBs, we assumed them 
to have equal luminosities. The results from 1000 sim¬ 
ulations are shown in Figure [8l where we plot the XLF 
using the mean (open circles; the standard deviation is 
shown as the error bar) of LMXBs in each luminosity 
bin from these simulations. The simulated XLFs seem 
steeper than the observed one overall, as expected. How¬ 
ever, the difference is small, which is mainly due to two 
reasons. One is that about one third of GC LMXBs are 
observed to cluster within a very narrow luminosity range 
(0.8-1.4)xl0 37 erg s _1 , but these GCs are widely spread 
around in the color-magnitude diagram, thus producing 
no significant source blending effects. The second reason 
is that in the area of parameter space where GCs prefer¬ 
entially host LMXBs, both faint and luminous LMXBs 
were observed. That is, the source blending effects are 
seen at different luminosity levels. Therefore we conclude 
that the source blending effects on the XLF are negligi¬ 
ble, within the uncertainties of our data. 

Instead of assuming equal weights to construct the 
XLF (Section 13. 4|) . we explored the option of weight¬ 
ing the incompleteness function to reflect the tendency 
of LMXBs to be detected in metal-rich and massive GCs. 
We recalculated the XLF of GC LMXBs with the incom¬ 
pleteness function weighted by the probability of hosting 
one or more LMXBs, based on Equation [3] again with 
A = 0.25 if M z < -10.0 and (g - z) > 1.13 and A = 0.14 
elsewhere (the source blending effects cannot be taken 
into account simultaneously though). We find that the 
XLF obtained in this way shows no significant difference 
from that shown in Section 13.41 

Therefore the observed paucity of faint GC LMXBs 
below the XLF break ~10 37 erg s _1 in NGC 3115 should 
be real. It could be explained if there is a transition 
from persistent sources to transi ents around this break 
HVosfi et. al .1120091:IKim et alJ200flh One main class of GC 
LMXBs could be ultracompact X-ray binaries (UCXBs). 
These are NSs accreting from white dwarf (WD) com¬ 
panions with very short orbital periods (P or b 1 hr) 



hi _L_ . I.l _i. .l 

1 0 36 1 0 37 1 0 38 

L x (erg s" 1 ) 


Fig. 8. — The XLF of observed GC LMXBs (filled black circles) 
and the XLF using the mean of 1000 simulations (open red circles, 
with the standard deviation as the error bar). The simulations 
take into account the possibility that some GCs might host mul¬ 
tiple LMXBs (see the text). The XLFs were not corrected for the 
incompleteness effects. 

and might be effectively produced in GCs through direct 
stella r c ollisions between NSs and red giants (jVerbuntl 
Il987jl . iBildst.en fc Delovel (|2004D first suggested that 
UCXBs with P or b ~ 8-10 min could explain the XLF 
of GC LMXBs at high luminosities. The second slope 
«2 = 1.6 ± 0.2 that we obtained in the broken PL fit to 
the XLF is consistent with their p r edicti on (02=1.77). 
However, according to lLasota et al.l (2 0081 ). for the tran¬ 
sition from persistent to transient behavior for a He-rich 
X-ray irradiated accretion disk to occur at around 10 37 
erg s _1 , systems with P or b > 40 min are preferred, at 
least below the break luminosity. 

It has been long-debated whether the entire popu¬ 
lation of LMXBs in galaxies, including those in the 
field, was formed dynamically in GCs dWh ite et al.1 
2002 IKundu et all [20021 [20071: fw(] 120051 l.luettl 120051 
Humphrey & Buote 120081) . Our result that the XLFs of 
GC and field LMXBs appear to be different agrees with 
previous findings fe.g. JVoss et al.l [20091 IKim et ahll2009l 
I Zhang et al .11201 H i . The difference indicates that they are 
formed through different channels that result in different 
system configurations (orbital period, mass ratio, etc.) 
and thus with different mass accretion rates. Therefore, 
our result supports the idea that not all field LMXBs are 
formed dynamically in GCs. 

In summary, source blending should occur in GC 
LMXBs, but we do not expect it to significantly affect 
the XLF. The observed paucity of faint GC LMXBs be¬ 
low the XLF break ~10 37 erg s _1 , compared with held 
LMXBs, is likely real, and one explanation is that GC 
LMXBs are dominated by accreting neutron stars with 
white dwarf donors that show a transition from persis¬ 
tent sources at high luminosity to transients at low lu¬ 
minosity around this break. The different XLFs of GC 
and held LMXBs suggest that held LMXBs are not all 
formed dynamically in GCs. 

5. CONCLUSIONS 
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We have studied the XLF of LMXBs in the early-type 
galaxy NGC 3115, using the the Megasecond Chandra 
XVP Observation of this galaxy. Including three previ¬ 
ous observations, we obtained a total exposure of ~ 1.1 
Ms and reached a detection sensitivity of L lim ~ 10 36 erg 
s _1 , much lower than typically achieved for other early- 
type galaxies (Lu m ~ 10 37 erg s _1 ) by Chandra. Our fit 
to the XLF of all LMXBs supports the presence of the 
low-luminosity break of the XLF at around 5.5 x 10 37 
erg s _1 , with the differential PL slope of 1.0 ± 0.1 and 
2.2±0.4 below and above the break, respectively, as seen 
in many previous stud i es (e. g, Gilfanov 2 0041 IVoss et al.1 
120091 : iKim et all 120091 : iZhang et all 201211 . However, we 
cannot exclude the possibility that the break is at around 
1.6 x 10 38 erg s _1 and is sharp, with the differential PL 
slope of 1.1 ± 0.1 and > 6.7 below and above the break, 
respectively. 

We further created the XLFs for field and GC LMXBs 
separately. Due to relative few GC LMXBs, the XLF 
of the field LMXBs is very similar to the XLF of all 
(GC±field) LMXBs and still shows degeneracy in the 
broken PL fit (one fit with a break at around 5.7 x 10 3 ' 
erg s -1 and the other fit with a sharp break at around 
2.0 x 10 38 erg s^ 1 ). The field LMXB XLF seems to show 
spatial variation, with the slopes and the break in the in¬ 
ner region ((0.046-0.2)L>2s) being smaller than those in 
the outer region (( 0 . 2 - 1 . 0 )-D 2 s)- This could be due to the 
incompleteness effects of the optical GC detection in the 
inner region and/or contamination of remnant LMXBs 
left behind from the destruction of GCs that drift to¬ 
ward the galactic center due to mass segregation. The 
XLF from the outer region is thus probably more close 
to the XLF of primodial field LMXBs. It has a differen¬ 
tial PL slope up to a break at around 1.7 x 10 38 erg s _1 , 
which is close to the Eddington limit of NS LMXBs. The 


detection of spatial variation explains the degeneracy in 
our fit to the XLF from the whole study region. 

The XLF of GC LMXBs overall is flatter than that of 
field LMXBs. Our observation of the difference between 
the XLFs of GC and field LMXBs casts doubt on the idea 
that all LMXBs in the galaxy are formed dynamically in 
GCs. The break of the GC LMXB XLF is at around 
1.1 x 10 37 erg s -1 and might be due to a transition from 
persistent sources at high luminosity to transients at low 
luminosity, which can be explained if GC LMXBs are 
dominated by accreting NSs with WD donors. 

As in previous studies, we found that metal-rich/red 
GCs are more likely to host LMXBs than the metal- 
poor/blue ones, an effect that is more significant for more 
luminous LMXBs, and that more massive GCs are more 
likely to host LMXBs. Although source blending is likely 
to occur, our simulations indicate that it should not sig¬ 
nificantly affect the XLF. 

Whi le at th e en d of the preparation of this paper, 
iLehmer et al.1 (I2014T ) also reported the study of three old 
normal galaxies. Their main goal was to test the evo¬ 
lutionary model of LMXBs, but they also obtained the 
XLFs of LMXBs in NGC 3115 using the same data pre¬ 
sented here. In the Appendix, we briefly compared our 
study with theirs. We found no large discrepancy be¬ 
tween our results and theirs, if factors such as the fitting 
degeneracy and the possible spatial variation are taken 
into account. 

The work is supported by Chandra XVP grant G02- 
13104X. This material is based upon work supported in 
part by the National Science Foundation under Grants 
AST-1211995 and AST-1308124. This material is based 
upon work supported in part by HST-GO-12759.02-A 
and HST-GO-12759.12-A. GRS acknowledges support 
from an NSERC Discovery Grant. 


APPENDIX 

COMPARISON WITH LEHMER ET AL. 2014 

In their study of the evolutionary model of LMXBs, ILehmer et al.1 (2014) also presented the XLFs of LMXBs in 
three normal galaxies, including NGC 3115. They used the same data presented here for NGC 3115, but there are 
many differences between their analysis method and ours. Lelnner et al. used the limiting significance level of 10 -5 
for wavdetect and kept sources with false binomial probability less than 0.004 (see their Equation 1), while we adopted 
the limiting significance level of 10 -6 for wavdetect (Paper I) and used sources above Ly lm , which results in using 
only sources with the signal to noise ratio >2.9. Moreover, Lelnner et al. used sources detected from the merged and 
individual observations (though the false binomial probability was calculated exclusively from the merged photometry), 
while we used sources detected from the merged observation only. Therefore we expect that Lehmer et al. could detect 
more real faint sources but also more spurious faint sources than we did. 

The region studied is also different. Lehmer et al. used an elliptical region of a semi-major axis of 2.7', a semi-minor 
axis of 1.1' and a position angle of 45° (based on the RT-band galaxy emission), excluding a central circular region of 
radius 20", but we use the slightly larger D 2 5 ellipse, excluding the central a = 10" elliptical region. We excluded a 
smaller central region because of our use of subpixel binning images for source detection in the central region. Within 
their study region, we have 90 sources above Lb m , among which 17.4 are expected to be CXB sources, while in our 
study region, we have 145 sources above Li im , 26.5 of which are expected to be CXB sources. The handling of the CXB 
contribution in the fits to the XLFs is also different. Lehmer et al. excluded all AGNs that they could identify from 
the LIST/ACS imaging (they fou nd 9 su ch sources , which is about 50% of the expected number), while we estimated 
the CXB contribution following iGeorgakakis et ahl (12008th 

We tried to check whether we can reproduce their XLF fitting results (they also fitted the XLFs with a broken PL) 
based on the sources that we detected but using their study region. We followed their technique to exclude all AGNs 
that we could i den tify (10 f rom Paper I) from the HST /ACS imaging instead of estimating the CXB contribution 
following IGeorgakakis et all ( 2008 ). We found that we generally obtained slightly lower values of the initial slope a.\ 
by (l-2)cr, most probably due to their inclusion of more very faint sources than were in our sample. Specifically, in the 
fit to the total XLF, they obtained a\ = 1.5 ± 0.1, and we obtained a.\ = 1.3 ± 0.1 if we chose a high-break solution 
with Lb = 1 - 571 q 53 x 10 38 er § S_1 similar to their Lb = (1.76 ± 0.02) x 10 38 erg s _1 . We still see the degenerate 
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low-break (Lb = (5.0 ± 0.2) x 10 37 erg s _1 ) fit, which has a C statistic larger than that of the high-break fit by only 
0.4 (or by only 0.1 if the CXB contribution was subtracted from modeling instead). For the field LMXB XLF, we 
preferred a high-break fit, which was also adopted by Lehrner et al. and is similar to our XLF of field LMXBs in 
the outer part of the D 25 region (Section [331). For the GC LMXB XLF, they also obtained a high break luminosity 
Lb = (1.76± 0.19) x 10 38 erg s _1 , but we cannot constrain it well (Lb = 4 . 5 + 3 4 g 2 x 10 37 erg s _1 ; we have 22 GC LMXBs 
(including three candidates) above Li; m in their region, while they have 25 GC LMXBs detected). Overall, we see no 
large discrepancy between our results and theirs for LMXBs in their study region. No large discrepancy is seen either 
between our results for our study region and theirs, if factors such as the fitting degeneracy and the possible spatial 
variation are taken into account. 
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